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t:j- ' Abstract 



A Bianchi type- VI cosmological model in the presence of a mag- 



^^ ' netic flux together with a cloud of cosmic strings is considered. In 

general, the presence of a magnetic field imposes severe restrictions re- 
garding the consistency of the field equations. These difficulties could 

K^ I be overtaken working either in a Bianchi type-VIo spacetime or assum- 

^ • ing a particular coordinate-dependence of the magnetic field. Using a 



few plausible assumptions regarding the parametrization of the cosmic 
strings, some exact analytical solutions are presented. Their asymp- 
totic behavior for large time is exhibited. 
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1 Introduction 

Since the observation of small anisotropics in the microwave background 
radiation (CMB) [1] and large scale structures [2] it became clear that a 
pure Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology could not 
explain all the properties of our Universe. It is therefore natural to consider 
anisotropic cosmological models that allow FLRW Universes as special cases. 

It is usually assumed that at the very early stages of the evolution of 
the Universe, during the phase transitions, the symmetry of the Universe 
was broken spontaneously [H H] . Topological defects such as strings, domain 
walls, monopoles has received considerable attention in cosmology since they 
could play an important role in the formation of large structure of the Uni- 
verse. 

String cosmological models have been used in attempts to describe the 
early Universe and to investigate anisotropic dark energy component includ- 
ing a coupling between dark energy and a perfect fluid (dark matter) [5]. 
Cosmic strings are one dimensional topological defects associated with spon- 
taneous symmetry breaking in gauge theories. Their presence in the early 
Universe can be justified in the frame of grand unified theories (GUT). 

A large number of astrophysical observations proves the existence of mag- 
netic fields in galaxies. Galactic magnetic fields which we observe today could 
be relics of a coherent magnetic field existing in the early Universe, before 
galaxy formation. Any theoretical study of cosmological models which con- 
tain a magnetic field must take into account that the corresponding Universes 
are necessarily anisotropic. Among the anisotropic spacetimes, Bianchi type- 
VI (BVI) space seems to be one of the most convenient for testing different 
cosmological models. 

The object of this paper is to investigate a BVI string cosmological model 
in the presence of a magnetic field due to an electric current together with 
the strings. Our paper is organized as follows: In Section 2 we derive the 
field equations of BVI cosmic string model in the presence of a magnetic field. 
Section 3 deals with the exact solutions obtained using some simple plausi- 
ble assumptions and describe their asymptotic behavior. The last Section 
contains conclusions. In the Appendix A, the geometrical properties and the 
shear tensor for BVI model are briefly described. 



2 Model and field equations 

The gravitational field in our case is given by a BVI metric: 

ds^ = dt^ - ale-^""' dx'^ - a^e^"^ dy^ - aj dz^ , (1) 

with ai, 02, as being the functions of time only. Here m, n are some arbi- 
trary constants and the velocity of light is taken to be unity. It should be 
emphasized that the BVI metric models a Universe that is anisotropic and 
space-dependent. The geometrical properties of the BVI spacetime are sketch 
in Appendix A, including the relationship with other Bianchi-type Universes. 
The Einstein field equations for BVI metric ([T]) are written in the form: 

^ + ^ + ^^-^ = .T,\ (2a) 

a2 (33 ^2 as «3 

as cii dgdi m^ 

— + — + 2 = ^^2 ^ (2b) 

as ax as ai a^ 

^ + ^ + ^^ + ^ = nTi, (2c) 

CLl 0,2 Oi\ 0,2 ao 

dido do dq ds di m^ — mn + n^ ^n , .s 

— — + —— + —— 2 = kT°, 2d 

Oi 02 02 as as ai ag 

<^l 0-2 / nOs „o /o \ 

m n [m — n) — = K-/3 . (2e) 

ai a2 as 

Here overdots denote differentiation with respect to time (t). The energy- 
momentum tensor for a system of cosmic strings and magnetic field is chosen 
to be 

T^ = pu^.u'' - Xx^x" + E"^ , (3) 

where p is the rest energy density of strings with massive particles attached 
to them and can be expressed as p = Pp + A, where Pp is the rest energy 
density of the particles attached to the strings and A is the tension density 
of the system of strings O O |8] which may be positive or negative. Here Ui 
is the four velocity and Xi is the direction of the string, obeying the relations 

UiU" = -XiX"" = 1, MjX* = . (4) 

For the electromagnetic field E^^ we adopt the form given by Lichnerowich 

m 

E^^ = fx[\h\'{uX-l5;)-h,h'']. (5) 



Here /x is the magnetic permeability and /i^ is the magnetic flux vector defined 
by 

h^, = -* F^^u", (6) 

H 

where *F^i, is the dual of the electromagnetic field tensor F^,^. 

In what follows the co moving coordinates are taken to he u^ = 1, u^ = 
u^ = u^ = 0. We choose the incident magnetic field to be in the direction of 
z-axis so that the magnetic flux vector has only one nontrivial component, 
namely h^ ^ 0. In view of the aforementioned assumption from ([6]) one 
obtains F23 = F31 = 0. We also assume that the conductivity of the magnetic 
fluid is infinite which leads to Fqi = -F02 = -fos = 0. Therefore there is only 
one non- vanishing component of F^^, namely F12. Then from the first set of 
Maxwell equation 

Fi_iu;p + Fyp-^^ + Fpfj,.^^ = , (7) 

where the semicolon stands for covariant derivative, one finds 



Then from ([6]) we get 



F12 = X, X = const. (8) 

hy, = exp Urn — n)z] . (9) 

Finally, for E^ one finds the following nontrivial components 

^0 = -El = -El = Ei = E= ^^ exp [2(m - n)z] . (10) 

Using comoving coordinates we have the following components of energy 
momentum tensor jllj : 

T° - p = -T,' = -Tl = Tl-\ = ^^exp [2(m - n)z] . (11) 

Taking into account that T3 = from ( l2el) one immediately finds 

r-^Y=xr-^y, at = const. (12) 

Let us now introduce a new function 

V = 010203. (13) 



Then from (1121) one finds 

a, = (^Ari;"ag^-2"j , (14) 

as = (t;"ar""/Ar) • (15) 

Summation of (12al) . pbj) . f l2cl) and three times fl2d|) gives 

- = 2 ^^ — + 3p + A + ^exp2m-nz). 16 

V ag 2 V yUf ^ / 

Let us note that from the energy-momentum conservation law one finds 

p+%-^X = 0. (17) 

V as 

Taking into account the ^-dependence of the energy-momentum tensor 
(fTTjl the r. h. s. of eqs. (j2]) have also a z-dependence, while the metric 
functions a^ depend on time only. Therefore, in general, eqs. (E]) are incon- 
sistent for a BVI model in the presence of a magnetic field. There are two 
possibilities to restore the consistency of eqs. (E]): 

1. to limit ourselves to Bianchi type-VIo (BVIq), namely, to consider the 
case 171 = 11] 

2. to assume a special ^-dependence of the magnetic permeability /i in 
order to compensate for the 2;-dependence of the magnetic flux com- 
ponent /is, eq. (jH]). A similar assumption was used by Bali [TU] in a 
different context 

In what follows we shall analyze these two possibilities in turn. 

3 Solutions of field equations 

3.1 BVIo model 

Assuming m = n we have 

ai = Ma2, M = const. (18) 



which permit us to express Oi and a2 in terms of 03 and v. 



a2 



M-^/\^/^a-"\ (20) 



Consequently Einstein's equations (|2]) are reduced to the following set of 
independent equations: 

(21a) 

(21b) 

(21c) 



03 'di (23 di m^ 
03 tti 03 ai a| 
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di d? m2 

~'~ 2 "*" 2 ~ 
tti a| ag 


— + A, 


d^ di ds m? 

2 ~'~ 2 ~ 
af ai 03 03 




where we introduced the notation 




^_.XW= 





2/i 



(22) 



Therefore there are three equations (1211) for four unknown functions Oi, 
03, p, A. It is customary to assume a relation between p and A in accordance 
with the equations of state for strings. The simplest one is a proportionality 
relation [6]: 

p = a\. (23) 

The most usual choices of the constant a are [TT], [T21 [131 [S] 

{1 geometric string 

1 + a; w > 0, p string or Takabayasi string (24) 

— 1 Reddy string. 

Using relation fl23|) between p and A we get that 

pva-^ = C, (25) 

with C an arbitrary constant. 

Now the system of differential equations is determined and we can proceed 
to solve it. However this system of nonlinear differential equations is quite 



intricate and we should resort to numerical simulations which will be reported 
elsewhere [15] . 

Here we limit ourselves to investigate the asymptotic behavior of the 
solutions for large t. For example, we can investigate the possibility to reach 
an isotropic regime, i.e. all functions Oj to have a similar behavior for t — > cxd 

ai ~ a2 ~ as ~ f 3 , (26) 

and consequently for the density of strings 



As it can be observed from (IA.11|) - flA.13|) in the case of isotropic space- 



time the components of the shear tensor a* vanish. 

To explore this possibility of isotropization, we shall investigate the equa- 
tion of evolution of v, ( TT6|) for m = n: 



V ^m^ K /„ , X^a3 



Assuming the asymptotic relation between Oj and v fl26l) we get from (1281) 
the following differential equation valid in the asymptotic regime t — > oo: 

V = Civ^ + C2V^ + C^v'^ , (29) 

allows the following first integral 

dv 



where Ci = 2m^,C2 = ^ 2a ' ^3 ~ ^ ^^^ some constants. This equation 



/ 



Civl + C^v^^^ + Cqv^ + C-r 



t + to , (30) 



where C4 = 3Ci/2, C5 = 6aC2/(3« + 1) and Cg = 3C3. Here to and C7 
are constants of integrations. 
We observe that for 

3^^ 3- (^^ 



i.e, a > 1 or a < , the term with f 3 is dominant in the integration ( 130|) 
and finally we get 

w ~ t^ (32) 



and consequently 

as ~ t , (33) 

and 

p~t-"^ (34) 

On the other hand for 

^>l. (35) 

i.e. a G (0, 1), the term with v^'^^^ is dominant in the integration ( 130|) and 
we obtain 

V ~ t3^ . (36) 

For a G (|, 1) we have a power growing in time for v and in the hmiting 
case a = I we get an exponential behavior in time. Finally, we note that for 
a e (0, |) there are no solutions in this model presenting an expansion of the 
Universe for large t. 

3.2 BVI model with a specific magnetic permeability 

As a second possibility to assure the compatibility of eqs. (E]) we assume a 
special z-dependence of the magnetic permeability [lO] : 



/i = /io exp [2(m — n)z] , (37) 

with fiQ a constant. Let us note that for z = the exponential factor is 1, 
but for z — y ±oo this factor vanishes or diverges, depending of the sign of 
the difference m — n. This unusual behavior of the magnetic permeability is 
accepted here as a working hypothesis. 

In this case eqs. ([2]) are compatible and could determine all unknown a^, p 
and A. As in the previous case, the numerical simulations will be presented 
elsewhere [15] and here we shall analyze only the asymptotic behavior of the 
solutions. 

Let us assume that for large t 



and consequently 





as 




ai 


rsj 


I) m+n ^ 
nj + m(l~2'Y) 
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r\^ 


1) m+n ^ 



(38) 



(39) 



and it is quite simple to verify that all equations ([2]) support this behavior 
in the asymptotic regime t — > cxd. For a particular value of 7, namely for 
7 = |, we recover the isotropization f l26|) discussed above. 

Let us observe that from eq. flTTl) we have in the asymptotic regime 



p = ^(-p + 7A). (40) 

We could consider that p has an asymptotic behavior correlated with that 
of V 

pr^v', (41) 

with 6 some constant which imply a proportionality relation between p and 
A as in eq. (^^, namely, A = ^^p. 

With all these assumptions a similar equation of evolution for v as in (130!) 
holds with appropriate constants of integration. The corresponding analysis 
of the asymptotic behavior of solutions of this equation proceeds as above. 

4 Conclusion 

We have studied the evolution of an anisotropic Universe given by a BVI 
cosmological model in presence of a cloud of cosmic strings and magnetic 
flux. It is found that the system with ^-dependent magnetic field within the 
scope of BVI spacetime is not consistent in general. This difficulty could be 
overcome working either in a BVIq metric, setting m = n, or introducing a 
particular z-dependence of the magnetic permeability. 

In a forthcoming paper [15] we shall present some numerical simulations 
and a detailed analysis of the stability and singularities of the field equations 
for the present cosmological model. 

Acknowledgments 

We are thankful to Suresh Kumar for helpful comments. This work is sup- 
ported in part by a joint Romanian-LIT, JINR, Dubna Research Project, 
theme no. 05-6-1060-2005/2010. M.V. is partially supported by GNOSIS 
program IDEI-571/2008 and NUCLEU program PN-09370102, Romania. 



Appendix A. BVI cosmological model 

The gravitational field in the present paper is given by a BVI cosmological 
model in the form ([1]). A suitable choice of m, n as well as the metric 
functions ai, 02, 03 in the BVI given by ([T]) evokes the following Bianchi- 
type Universes: 

• for m = n the BVI metric transforms to a BVIq one, i.e., m = n, 
BVI =^ BVIq G open FRW with the line elements 

ds^ = df - ale-^""' dx^ - a^e^"^ dy'^ - a^ dz^ ; (A.l) 

• for m = —n one gets the Bianchi-type V (BV) spacetime; 

• for n = the BVI metric transforms to a Bianchi-type III (Bill) one, 
i.e., n = 0, BVI =^ Bill with the line elements 

ds^ = d^ - ale-^"^^ dx^ - al dy^ - 4 dz^ ; (A.2) 

• for m = n = the BVI metric transforms to a Bianchi-type I (BI) one, 
0, BVI =^ BI with the line elements 

ds'^ = dt^ — a\ dx^ — a\ dy'^ — a^ dz^ ; (A. 3) 

• for m = n = and equal scale factor in all three directions the BVI 
metric transforms to a Friedmann-Robertson- Walker (FRW) Universe, 
i.e., m = n = and a = h = c, BVI =^ FRW with the line elements 

ds'^ = dt^ - c? [dx^ + dy^ + dz^) . (A.4) 

Let us go back to the BVI cosmological model ([1]). The nontrivial Chris- 
toffel symbols of the BVI metric read 

2 2 

pi _ _™ F^ - n F^ - !^p-2m2 p3 _ ^^2 2nz 

^ ZX ~ '"'' -*- 32 ~ '''5 -*- 11 ~ 2 ' ■*■ 22 ~ 2 

ag ag 
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The non- vanishing components of Riemann tensor corresponding to ([T]) 
are 



a2 as 
02 as ' 
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ai\ 
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The nontrivial components of the Ricci tensor are 

dO / ai 02 . -as\ 

Kn = —im n [m — n) — 1 

V tti 02 03/ 

Voi 02 as 

1 f'di (21(22 (2i (2s m? — mn 
R^ = \ \ 2 

\ai Oi 02 ai as ag 

2 /a2 (2i (22 (22 (2s 71^—17171 

R2 = I 1 1 n 

Va2 ai 02 a2 as ag 

3 /as (2i(2s , (22 (2s ITl'^ + n"^ 

rin = ( 1 1 K 

Vas ai as a2 as a^ 

The Ricci scalar reads 

D_2r^^ ^^ "^ '^^ '^^ "'^^^ ^^^^ m^-mn + n^ . 
ai a2 as ai a2 a2 as asai a| 

Let us now find expansion and shear for BVI metric. The expansion is 
given by 

^ = <, = < + r;:„«° , (A.6) 

and the shear is given by 

a' = ^a,,a^\ (A.7) 

with 

a^u = \ W;aP: + M.;aP;] " \^P^.u , (A.8) 

11 



where the projection vector P: 

p2 = p, p^^ = g^^ - u^u,, pi: = 6^- U^U, . (A.9) 

In comoving system we have -u^ = (1, 0, 0, 0). In this case one finds 



^ = ^ + ^ + ^, (A.IO) 

(Xl 0,2 fl3 



and 



„„ = !*^(_2:^ + :k + :iL) ^ ,j . 4 (_2ii + :k + :^) , (A.u) 

6 V ai a2 03/ o V ai a2 03/ 

,,, = Mf!:(_2:kH.:^ + M^,._ i(_2:^ + :k + ii), (A.12) 

o V 02 03 Oi / o V 02 a3 Oi / 

ag / ^03 ^ di ^ 02^ 3 1 / ^^3 , di , 02 



^33 = :^ -2::^ + - + -^a3^ = -- -2::^ + - + - . (A.13) 
o V 03 ai 02/ oV 03 fli 02/ 
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